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Abstract 



Using the formalism of relativistic acoustic geometry we study the expanding chiral 
fluid in the regime of broken chiral symmetry near the QCD chiral phase transition 
^■i temperature T c . The dynamics of pions below T c is described by the equation of motion 

(— 1 \ for a massless scalar field propagating in curved space-time similar to an open FRW 

universe. The metric tensor depends locally on the soft pion dispersion relation and the 
. four- velocity of the fluid. In the neighborhood of the critical point an analogue trapped 

region forms with the analogue trapped horizon as its boundary. We show that the 
associated surface gravity diverges near the critical point as k ~ (T c — T)" 1 . Hence, if 
OO 1 the horizon forms close to the critical temperature the analogue Hawking temperature 

may be comparable with or even larger than the background fluid temperature. 

o; 

1 Introduction 

Analog gravity models of general relativity seem promising routes to providing laboratory 
tests of the foundation of quantum field theory in curved space-time [1] (for a review and an 
extensive list of references, see [2]). The models are based on kinematics of waves propagating 
in inhomogeneously flowing media. Various aspects of these phenomena have been studied 
in acoustics [3] optics [4] and superfluidity [5]. In this paper we study the analogue gravity 
model based on massless pions propagating in a hadronic fluid. 

Strongly interacting matter is described at the fundamental level by a nonabelian gauge 
theory called quantum chromodynamics (QCD). At large distances or small momenta, the 
QCD exhibits the phenomena of quark confinement and chiral symmetry breaking. Chiral 
symmetry breaking expresses the fact that massles quarks confined in hadrons appear effec- 
tively as massive constituents with a dynamically generated mass of several hundred MeV. 
At low energies, the QCD vacuum is characterized by a nonvanishing expectation value [6, 7]: 
(ipip) ~ (235 MeV) 3 , the so called quark condensate, which describes the density of quark- 
antiquark pairs found in the QCD vacuum and its nonvanishing value is the manifestation 
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of chiral symmetry breaking. The chiral symmetry is restored at finite temperature above 
certain critical temperature T c of the order of 150 MeV. In the chirally broken phase, below 
T c the pions, although being massless, propagate slower than light [8, 9, 10] with velocity 
approaching zero at the critical temperature. The pions propagate through the expanding 
fluid the velocity of which may well be comparable with the speed of light. Hence, it is very 
likely that the flow velocity will exceed the pion velocity if the temperature of the fluid is 
below and close to the critical temperature. A similar condition is also realized in cosmolog- 
ical models of expanding universe in which the recession velocity of the expanding spacetime 
beyond the Hubble horizon exceeds the speed of light. 

In order to explore this analogy with cosmology more closely we need to specify the 
velocity field of the flow in a given spacetime geometry and the fluid parameters such as the 
density, temperature and the velocity of pions propagating in the fluid. For that purpose we 
exploit a simple boost invariant Bjorken type spherical expansion [11] of the chiral fluid. A 
similar model has been studied some time ago in the context of disoriented chiral condensate 
[12]. Furthermore, the velocity of pions in the chiral fluid can be derived using a linear 
sigma model as an effective low energy model of strong interactions. Then, the propagation 
of massless pions provides a setting for a geometric analogue of expanding space time. Since 
the velocity of the fluid may exceed the velocity of pions it is conceivable that an apparent 
horizon may exist with the associated analogue Hawking radiation. In this paper we study 
the analogue Hawking effect using the Kodama-Hayward definition of surface gravity [13]. 

In the next two sections we describe our basic concepts and in sec. 4 we derive an 
expression for surface gravity and study its behaviour near the critical point. 



2 Dynamics of the chiral fluid 

In the following we assume that the chiral fluid undergoes a boost invariant Bjorken type 
spherical expansion [11]. In this model the radial three- velocity in radial coordinates x M = 
(t, r, <f) is a simple function v = r/t. Then the four- velocity is given by 

u " = (t/r,r/T,0,0), (1) 

where r = \/t 2 — r 2 is the proper time. With the substitution 

t = r cosh y, 

r = rsmhy, (2) 



v = tanhy, (3) 



the radial velocity is expressed as 
and the four-velocity as 

if = (cosh y, sinh y, 0, 0) . (4) 

The substitution (2) may be regarded as a coordinate transformation from ordinary radial 
coordinates to new coordinates (r, y, ip) in which the flat background metric takes the form 

= diag (l, -r 2 , -r 2 sinh 2 ?/, -r 2 sinh 2 y sin 2 6>) , (5) 
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and the velocity componets become = (1,0,0,0). Hence, the new coordinate frame 
is comoving. The metric (5) corresponds to an FRW expanding cosmological model with 
cosmological scale a = r and negative spatial curvature. 

The temperature of the expanding chiral fluid, to a good approximation, is proportional 
to t -1 . This follows from the fact that the chiral matter is dominated by massless pions, and 
hence, the density of the fluid may be approximated by the density p = (gir 2 /30)T 4 of an 
ideal massless boson gas [14]. Using this and the energy-momentum conservation one finds 

T=^, (6) 

T 

where the constant Co may, in principle, be fixed from the phenomenology of high energy 
collisions. 

Next we derive the pion velocity in the chiral fluid. The dynamics of mesons in a medium 
is described by a chirally symmetric Lagrangian of the form [9, 10, 15] 

1 Tft 2 \ 

C = -{agT + buV)d^d uV - " - A ^ 2 )\ (7) 

where is the velocity of the fluid, and g^ is the background metric. The mesons tp = (a, 
7r) constitute the (|, |) representation of the chiral SU(2) xSU(2). The parameters a and b 
depend on the local temperature T and on the parameters of the model m and A and may 
be calculated in perturbation theory. At zero temperature the medium is absent in which 
case a = 1 and b = 0. 

If itLq < the chiral symmetry will be spontaneously broken. At the classical level, the 
a field develops a nonvanishing expectation value such that (a) — f n . At nonzero tempera- 
ture the expectation value (a), usually referred to as the chiral condensate, is temperature 
dependent and vanishes at the chiral transition point. Hence, the quantity (a) serves as an 
order parameter. For temperatures below the chiral transition point the meson masses are 
given by 

ml = 0; m 2 a = 2\(a) 2 , (8) 

in agreement with the Goldstone theorem. The temperature dependence of the chiral con- 
densate (a) is obtained by minimizing the thermodynamical potential Q = — (T/V)ln Z 
with respect to (a) at fixed temperature T. At one loop order the solution to the extremum 
condition as a function of temperature exhibits a weak first-order phase transition [16, 17]. 
However, Pisarski and Wilczek have shown on general grounds that the phase transition in 
SU(2)xSU(2) chiral models should be of second order [18]. Hence, it is generally believed 
that a first-order phase transition in this case is an artifact of the one loop approximation. 
Two loop calculations [19] make an improvement and confirm the general analysis of [18]. 
Propagation of pions is governed by the equation of motion 

■^=d li [^ = g(ag' u/ + bu' i u v )d v ir] + V(a, = 0, (9) 

which follows from the effective Lagrangian obtained from (7) by spontaneous chiral sym- 
metry breaking. The quantity V in (9) is the interaction potential the form of which is 
irrelevant for our consideration. From (9) we obtain the pion velocity squared as [15] 

^ - Tti (10) 
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The parameters a and b at nonzero temperature may be derived from the finite temperature 
self energy £(g, T) of the pion in the limit when the external momentum q approaches 0. 
For a flat background geometry g^ u = r]^, the inverse pion propagator A -1 is given by 

A" 1 = aq»q» + b(q^) 2 - m 2 , (11) 

or in comoving frame, i.e., in a reference frame in which uF 1 = (1, 0, 0, 0), 

A" 1 = (a + b)q 2 - aq 2 . (12) 

Hence, the parameters a and b, can be expressed in terms of second derivatives of S(g, T) 
with respect to q and ^ evaluated at q^ = 0. 

The pion velocity at nonzero temperature was calculated at one loop level by Pisarski 
and Tytgat in the low temperature approximation [8]. At one loop level the only diagram 
that gives a nontrivial g-dependence of E(g, T) is the bubble diagram and the calculation 
may be performed for the hole range of temperatures below the chiral critical point [10, 15]. 
As we are particularly interested in the behaviour near the critical point of the chiral phase 
transition we make use of the exact results near the critical temperature based on scaling 
and universality analysis of Son and Stephanov [9] . In their notation our parameters a and 
b are expressed as 

f 2 f 2 

a+6=-% (13) 



<<r>" (ay 

where f s is equal to the temperature dependent pion decay constant and f 2 coincides with 
the isospin susceptibility Xi5- The chiral condensate (a) up to a factor of the order f 2 equals 
the quark condensate (ipip) at zero quark masses. In the limit T — > T c the quantity f 2 goes 
to a nonzero constant whereas, contrary to naive expectations, the scaling of f s turns out to 
be different than that of the order parameter (a): 

/ S 2 ~(T C -T)^; (a)~( Tc -Tf, (14) 

where d is the number of space dimensions and v and j3 are positive critical exponents. 
Hence, we obtain 

a~ (T c -T)( d - 2 ^ 2/3 ; a + b ~ (T c - T)~ 2/3 , (15) 
so the pion velocity scales as 

c 2 ~(T c -T)( rf - 2 ^, (16) 

as one approaches the critical temperature. In d = 3 dimensions, the critical exponents for 
the 0(4) universality class are v = 0.73 and (3 = 0.38 [20]. 



3 Chiral geometry 

We now introduce the analogue gravity metric which describes the effective geometry of the 
expanding chiral fluid using the the formalism of relativistic acoustic geometry [21, 22, 23]. 
The equation of motion (9) may be written in the form 

-1= d,(V^G GT)d v n + ^V(a, k)tz = 0, (17) 
v— G a 
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with the analogue metric tensor, its inverse, and its determinant given by 



Gp,v = —[g^u - (1 - 4) t V i "]> ( 18 ) 

C77- 



sT - (i - 4)*^ 



c 2 



(19) 



G = ^. (20) 



Hence, the pion field propagates in a (3+l)-dimensional effective geometry described by the 
metric G^ u . It is convenient to work in comoving coordinates (r,y, $,</?) with background 
metric defined by (5). In these coordinates the analogue metric tensor (18) is diagonal 
with components 



Gp, = — diag (4, -r 2 , -r 2 sinh 2 y, -r 2 sinh 2 ?/ sin 2 fl) , (21) 



where the parameters a and c n are functions of the temperature T which in turn is a function 
of r by (6). In the following we assume that these functions are positive. 

In contrast to [22], where it was assumed that both the background geometry and the 
flow were stationary, in an expanding fluid the flow is essentially time dependent. Hence, 
the acoustic geometry formalism must be adapted to a non-stationary space time. To this 
end we introduce the concept of analogue marginally trapped surface or analogue apparent 
horizon. 

To define the analogue apparent horizon we need to examine the behaviour of radial null 
geodesies of the analogue metric (21) in which a and c n are functions of r. In the following 
we assume spherical symmetry and denote by 1+ and It the vectors tangent to outgoing and 
ingoing affinely parameterized radial null geodesies normal to a sphericall two-dimensional 
surface. The tangent vectors are null with respect to the metric (18), i.e., 

Gfxvl+l 1 ^. = G^VtM^ = 0. (22) 

Using the geodesic equation 

i^v^r = 0, (23) 

where the symbol denotes a covariant derivative associated with the metric (18), one 
easily finds the tangent null vectors corresponding to future directed radial null geodesies, 

/^ = -fl,±-,0,0V (24) 



ar \ t 



The null vectors /+ and It point towards increasing and decreasing y, respectively. Hence, 
we adopt the usual convention and refer to 1+ and It (and the corresponding null geodesies) 
as outgoing and ingoing. 

The key element in the study of trapped surfaces is the expansion parameter e± of null 
geodesies. A two-dimensional surface S with spherical topology is called a trapped surface if 
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the families of ingoing and outgoing null geodesies normal to the surface are both converging 
or both diverging. More precisely, the expansion parameters 



e± = V„Z£ (25) 

on a trapped surface S should satisfy > 0. A two-dimensional surface H is said to 

be future inner marginally trapped if the future directed null expansions on H satisfy the 
conditions: s+\h = 0, Itd^e^u > and S-\h < 0. We shall refer to this surface as the 
apparent horizon since it is equivalent to the apparent horizon in cosmo logical context. 
From (25) with (24) we find 

From (26) one finds the condition for the apparent horizon 

— ± — V > = 0. 27 

v 



\foT/a K 



This equation defines a hypersurface which we refer to as the analogue trapping horizon. 
Any solution to Eq. (27), e.g., in terms of r for an arbitrary fixed t — i*, gives the location 
of the analogue apparent horizon th- 

The critical behaviour of the null expansions e + and e_ and of the derivative l^Ld^e + near 
the trapping horizon is of particular interest. Using (14)- (16) in the neighborhood of the 
horizon we find 

£+ ^(£l_J^), (28 ) 
ar£ \v t — t c / 

- - 4 " (29) 



ar c (r - r c ) ' 

r_ 9 „ £+ , '"-y*^) , (30) 

a 2 T^(r-T c ) 2 

where 

ri = 0-(d-2)v/4. (31) 

and r c = c /T c is the critical proper time corresponding to the critical temperature T c . Note 
that for an arbitrary fixed t = t* the outgoing null expansion e + vanishes at the radius r = r fl- 
at which v = c n (r — t c )/(i]t c ). The ingoing null expansion is obviously negative whereas the 
derivative of e + given by (30) is positive at r H and, according to the standard covention [24], 
the region {r > rn,t — i*} is future trapped with the apparent horizon located at ru as its 
inner boundary. 
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4 Surface gravity and analogue Hawking effect 



On suggestion by Hajicek [25] it was recently argued [26, 27, 28, 29, 30, 31] that the Hawking 
effect might be associated with the apparent horizon rather than the event horizon. This 
observation is irrelevant for a stationary space time where the apparent and event horizons 
coincide. In this case the apparent horizon is Killing and the surface gravity is uniquely 
defined as a parameter that measures the inaflinity of the properly normalized Killing vector 
The surface gravity k of the Killing horizon can be defined by 

TV^ = (32) 

evaluated on the horizon. If the geometry were stationary, i.e., if the components of G^ v in 
(18) were be time independent, the apparent horizon would coincide with the chiral event 
horizon at the surface defined by v = c n . In that case the surface gravity would read [22] 

*=l^!>-0, (33) 

where the derivative is to be taken at the horizon. 

In the case of nonstationary spacetime, the apparent horizon is neither Killing nor null. 
The definition of surface gravity in this case is not unique [28] and several ideas have been 
put forward how to generalize the definition of surface gravity for the apparent horizon 
[13, 24, 26, 32]. In this paper we addopt the prescription of [13] which, we believe, is most 
suitable for spherical symmetry. This prescription involves the so called Kodama vector 
[33] which generalizes the concept of the time translation Killing vector to non-stationary 
spacetimes. The Kodama vector we define as [24, 34] 

K a = ke al3 n/3 for a = 0, 1; IC = for i = 2, 3, (34) 

where e a/3 is the covariant two-dimensional Levi-Civita tensor in the space normal to the 
surface of spherical symmetry and n a is a vector normal to that surface. For the metric (21) 
n a is given by 

n M = <9 M ^Ty^sinhyJ . (35) 

Our definition differs from the original one [24, 34] by a normalization factor k which we 
have introduced in order to meet the requirement that should coincide with the time 
translation Killing vector £ M for a stationary geometry. Using the definition (34) with metric 
(21) we find 



k 9 T (t 



K T = cosh K y = —k — -sinh?/, (36) 

x/aCir ra 



with norm squared 

\K\ 2 = h afi K a K^ = k 2 smh 2 y 



coth 2 y 



(37) 
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The vector K a is spacelike in the trapped region and vanishes on the trapping horizon. 

In order to fix k, consider a general spherically symmetric analogue geometry with metric 
(18). In (t,r) coordinates 

n M = <9 M (r^J a/c^j . (38) 

Using (34) with metric (18) we find 



K l = ka 1 d r (^\Ja~Jc^j ; K r = —ka 1 rd t \ja/c K . 



(39) 



In the case of stationary geometry K r = because the quantities a and c w depend only on 
r. Then, if we set 

k = y/ac^ ^1 + ^Jc^/ardr^Ja/c^ , (40) 
we have = ^ in the stationary case. In (r, y) coordinates 

/ 



cosh y — - smh y \ . 



V 



In analogy with (32) the surface gravity k is defined by [13, 35] 

K a V [a K p] = K,Kp, 

or equivalently by 



(41) 



(42) 



(43) 



2^ 

where the quantities on the righthand side should be evaluated on the trapping horizon. 
Using (38) and (40) it is easy to check that this expression reduces to (33) for a stationary 
spherical flow, i.e., a steady flow where the velocity v and the constants a, b, and c n are time 
independent. By making use of the horizon condition (27) for our time dependent Bjorken 
spherical flow we find from (43) 



' = \—^==d n (\/^-hh a Pknp), 



K 



yi^ 2 

(i + c w vy 



c w | cl-clv-c^v 2 civ 2 | d 2 T (r^a~/Z 



2tv 



2t 



'aCj, 



(44) 



where v = tanhy and it is understood that the righthand side should be evaluated on the 
trapping horizon. 

The critical behaviour of k is determined by the behaviour of the quantity \J a/c v in the 
neighborhood of the critical point. From (14)-(16) we have 



oc (r - re)"" 



(45) 



in the neighborhood of the critical point, where rj is defined in (31). Then, on the trapping 
horizon in the limit r — > r c the two dominant terms in square brackets in (44) scale as 



Cn_ <9r (rygg I 



(46) 



V V ' ' ~ (77 + l)(r - r c )-\ (47) 

Plugging these expressions in (44) we find that the surface gravity in the neighborhood of r c 
behaves as 

k~ (r ] + l/2)(r-r c )- 1 . (48) 

With the critical exponents v = 0.73 and f3 = 0.38 for the 0(4) universality class in d = 3 
dimensions [20] we find 77 = 0.1975. 
The temperature 

T H = £ (49) 

associated with the analogue apparent horizon is the analogue Hawking temperature of 
thermal pions emitted at the apparent horizon as measured by an observer near infinity. 
Since the background geometry is flat, this temperature equals the locally measured Hawking 
temperature at the horizon. From (48) it follows that the analogue Hawking temperature 
may be arbitrary large in the limit when the analogue horizon approaches the critical point. 
However, Eq. (48) holds for an ideal spherical expansion of the chiral fluid and in a realistic 
scenario one could only expect that the Hawking temperature may be comparable with or 
larger than the background fluid temperature. 



5 Discussion 

It is tempting to speculate about possible signals for the analogue Hawking effect in a 
hadronic fluid. In principle, one could measure the temperature by fitting the pion spectrum 
to the thermal Planck distribution. However, one must be cautious before making any pre- 
diction for a realistic physical system. First, one must invent a reliable signal to distinguish 
between the thermal pions produced above the critical temperature from those emitted as 
an analogue Hawking radiation from the apparent horizon below the critical temperature. 
Second, a spherically symmetric expansion model considered here is not realistic for high 
energy heavy ion collisions. A more realistic hydrodynamic model would involve a trans- 
verse expansion superimposed on a longitudinal boost invariant expansion. In this case the 
calculations become rather involved as the formalism for general nonspherical spacetimes is 
not yet fully developed. This work is in progress. 
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